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1. Introduction 

Description of the nucleon-nucleon interaction, the basic ingredient of nuclear physics, di- 
rectly from quantum chromodynamics (QCD) is a daunting task. The breakdown of the pertubative 
expansion in terms of the QCD gauge coupling g necessitates the use of alternative methods. One 
such method proposed by 't Hooft in 1974 is to consider a QCD-like theory with the number of 
colors N c and the gauge group SU (N c ), large N c QCD [|3j]. The observables in large N c QCD are 
expanded in powers of l/N c around the large N c limit, g — > 0, N c — > °° and finite g 2 N c . In addi- 
tion, it is assumed that large N c QCD is a confining theory and the asymptotic states are SU (N c ) 
singlets. Despite our inability at present to evaluate even the leading order terms, a great deal 
of insight comes from knowledge of the scaling of hadronic observables in powers of 1 /N c . The 
phenomenological implications of large N c QCD are essentially topological in nature. 

The description of the meson and baryon observables in the large Nc limit requires different 
methods. Formally, it is due to the fact that the correlation functions in the meson sector have a 
smooth expansion in powers of 1 /N c while the correlation functions in the baryon sector diverge in 
the large N c limit. Physically, it is due to the fact that as shown by 't Hooft, the QCD in the large 
N c limit is a theory of an infinite number of stable non-interacting mesons. 1 The meson masses and 
sizes are independent of N c . This picture is phenomenologically satisfactory since in the real world 
mesons interact weakly. 

QCD is also a theory of strongly interacting baryons, the states carrying quantum numbers of 
the odd number of quarks. As was shown by Witten a consistent large N c description of strongly 
interacting baryons is possible. Remarkably, the same feature which on the one hand makes an 
analysis of the baryons in the large N c limit far more challenging than that of mesons, on the other 
hand allows one to apply a well-known method of nuclear physics, the semiclassical mean-field 
theory. Indeed, as argued by Witten, baryons in the large N c limit contain N c quarks with n-quark 
force scaling as N^". Thus, one can treat this weakly interacting many-body state in a mean-field 
approximation. Unfortunately, the explicit treatment is only available for heavy non-relativistic 
quarks, in which case the mean-field treatment corresponds to the Hartree approximation The 
picture that arises from such a treatment is that of a baryon with a mass of order of N c and a size 
and shape which are independent of N c . Despite the fact that explicit mean-field treatment in the 
case of the light quarks is unknown, the large N c scaling for baryon observables is expected to be 
valid. 

Witten realized [Q] that the above scaling of mesons and baryons indicates that the baryons 
in large N c limit arise as quantized soliton-like configurations of mesonic fields. A particular 
model which satisfies the large N c scaling is a well-known Skyrme model [|||. In this model the 
baryons appear as quantized skyrmions, the topological solitons of a particular non-linear mesonic 
lagrangianj^]]. The stability of baryons as quantum solitons is due to the existence of conserved 
topological current []8j 

In addition to a single-baryon sector, it is of great interest to consider baryon-baryon interac- 
tion in the large N c limit. Since the baryon mass diverges in the large N c limit the baryon-baryon 
scattering observables don't have a smooth limit for scattering at fixed center-of-mass energy and 
momentum transfer. For such momentum, p ~ N c , one instead focuses on the potential between 

1 Additionally, the spectrum of large N c QCD contains glueballs with vanishing mixing to mesons 
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two baryons which scales as N c . In this context, the large N c scaling rules had been used to analyze 



the spin-flavor structure of the nucleon-nucleon potential fllOQ. In particular, it was shown that at 



leading order the nucleon-nucleon potential is symmetric under contracted SU (4) spin-flavor sym- 



metry Jlip. In addition, one can also address a question of consistency of the meson-exchange 
picture of the nucleon-nucleon potential Ji2| , Up. 

The focus of the present talk is on the kinematic regime corresponding to a fixed center-of- 
mass velocity. In this case, both the kinetic and potential energy of two baryons are of order N c , 
and thus one expects a smooth limit to exist for the scattering cross-section. As argued by Witten 
[Q], here, as in the case of the single-baryon sector, the appropriate framework is the mean-field 
description. However, in this case one has to use time dependent mean-field theory (TDMFT). 
As was shown in [Q], one can discuss the spin-flavor structure of the total nucleon-nucleon cross 
section at leading order in large N c expansion. As in the case of the nucleon-nucleon potential, the 
emergent contracted SU (4) spin-flavor symmetry leads to certain relations between total proton- 
proton and proton-neutron cross sections. It will be shown in section || that these relations satisfy 
the behavior of the experimental nucleon-nucleon cross sections at the center-of-mass energies of 
order of a few GeV. 

2. Time Dependent Mean Field Theory Framework 

The goal here is to show that time-dependent mean-field theory (TDMFT) framework valid in 
the large N c limit, and contracted SU (2Nf) symmetry, where Nf is the number of light quark fla- 
vors, to make model-independent predictions about the spin-isospin structure of the total nucleon- 
nucleon cross sections. As discussed above, TDMFT treatment is a valid framework for the 
nucleon-nucleon scattering when the center-of-mass transfer momentum is p ~ N c . Since the nu- 
cleon size and hence the size of the interaction region are of order of Af°, the scattering in this 
kinematic regime is semi-classical. 

The description of interaction by TDMFT methods requires time-averaging over all field con- 
figurations consistent with the initial state of two nucleons. This precludes one from being able 



to calculate the S-matrix elements [|14j, |15|]. However, as shown in [J2p , there are certain inclusive 
nucleon-nucleon observables which can be evaluated in TDMFT framework. One such observable 
can be formed from conserved baryon current whose expectation value in the initial two-nucleon 
state can be in principle evaluated. The expectation values of the baryon current can be related to 
the inelastic differential cross section. 

In TDMFT framework each quark and gluon field of two nucleons move in a time-dependent 
field created by all other quarks and gluons. These equations are not known explicitly. As a result, 
one can not determine the nucleon-nucleon cross section even in the large N c limit. However, 
it is possible to determine the spin-isospin dependence of the cross section using the contracted 
SU (4) symmetry valid in the large N c limit. Since the focus here on spin-isospin dependence of 
the nucleon-nucleon scattering, one can use the Skyrme model which encapsulates the spin-flavor 



structure of large N c baryons [16]. 

In the Skyrme model the nucleon dynamics is described in terms of classical soliton config- 
urations built out of pion fields. A convenient form for such a soliton is given by SU (2)-valued 
matrices Uh{r) = exp(iTnF(r)), where % a are Pauli matrices, h = r/r and F{r) is the magnitude 
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of the pion field, F(r) = | ?f (r) | ^J§. Such classical configurations impose correlations between 
spacial and isospin rotations and they are referred to as hedgehogs. The baryons appear after the 
quantization of classical hedgehogs. This is done by quantizing the slow rotation of the hedgehog 
in isospin space, A^(t)Uh(r)A(t), given by the time-dependent SU (2) matrix A(t). These rotations 



describe the slow collective degrees of freedom of the hedgehog, the zero-modes [ |17| , |18Q , After 
quantization, the generators of these vibrations proportional to Tr (^x a A{t)A^ (t) \ correspond to the 
spin and isospin quantum numbers of the the ground-state band of states, I = J = 1/2, 3/2, ...N c /2 
in the large N c limit. The first two states correspond to nucleon and A baryons. The masses of these 
states are degenerate up to the terms /(/ + 1 ) /Mb ~ 1 /N c . This is a representation of the contracted 
spin-flavor symmetry SU (4) in the context of the Skyrme model. The spin-isospin dependence 
of the wave-function of these states is given by Wigner matrices 0^,(A) where m,m' are the 
third components of spin and isospin respectively. Note that the Wigner matrices are functions of 
parameters of collective rotations A and not of A(t). 

A crucial consequence of the above semi-classical analysis of the single-baryon sector in the 
Skyrme model is the appearance of the scale separation in the dynamics of the collective degrees of 
freedom described by the Wigner matrices Oj~i(A) at leading order in l/N c and intrinsic degrees 
of freedom which describe all other non-collective excitations which include excited states of the 
ground states baryons and mission of virtual and real mesons. The frequency of the collective 
excitations are of order 1 /N c while that of the intrinsic excitations are of order N®. This scale 
separation enables an adiabatic or treatment of the collective degrees in the context of TDMFT 
treatment of the nucleon-nucleon scattering analogous to the Born-Openheimer approximation in 
the context of the rotational and vibrational excitations of molecules. There the slow degrees of 
freedom correspond to vibrations of atomic nuclei in the averaged field produced by electrons 
whose motion represent the intrinsic excitations. 

To obtain an observable describing the nucleon-nucleon scattering one can start with a function 
which describes the initial state of two well separated hedgehogs corresponding to the initial state 
of two nucleons. As discussed in [0], in the context of the Skyrme model it is convenient to 
choose a conserved baryon current (?,t;A\ ,A2, v,S,w^ . In the Skyrme model it is a topological 
invariant. The dependence of the current on the collective degrees of freedom are described by the 
variables, A\,A2,v,b,n where A\,A2,v,b,n define the spin and isospin configuration, the center-of- 
mass speed, impact parameter and the unit vector along the direction separating the centers of the 
two hedgehogs at the initial moment. The initial distance is not indicated. Such parametrization 
corresponds to the semi-classical description of the scattering which as discussed above is valid 
in the large N c limit. The functional dependence of this current can only be determined once the 
explicit form of TDMFT equations is known. However, as shown below one does not need to know 
these equations to determine the spin-isospin structure of the corresponding scattering observable. 
This structure is determined by transformational properties of the current under the spin and isospin 
rotations which is determined by the contracted sin-flavor symmetry. 

As shown in [JjJ ||], the classical current can be turned into a differential cross section by 
integrating the current J^ 1 ( r , t ; A i , A2 , v, b , h ) over time and the impact parameters, 



dOin C (v,Ai,A 2 ;d,<l>) _ 2 



dQ. 



limfl 2 [ db(2nb) [ dtr(6,<l>) ■j(Rr(a),t;A 1 ,A z ,v,b,n\ , (2.1) 
Jo Jo V / 
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The above equations gives the probability for one hedgehog to emerge in a cone with a solid angle 



dQ. around a direction given by polar angles 6 and (p. The current in Eq. ( |2.1| ) is normalized as to 



give the total baryon number two. The time t = in Eq. ( |2.1| ) corresponds to the time at which two 



hedgehogs have the smallest separation. The integral in Eq. (gJJ) can be explicitly evaluated only 
when TDMFT equations are known. It is also important to note that the probability in Eq. (2.1) is 
integrated over all outgoing meson degrees of freedom. In the final analysis it will give the inelastic 
cross section. 

To turn the probability in Eq. ( |2.1[ ) into a nucleon-nucleon cross section one needs to evaluate 
an expectation value of do mc /dQ. in an initial two-nucleon state described by spin and isospin 
projections m\,m\ and ra2,n4 on the direction given by the unit vector h. It can be done due 
to the scale separation between the collective and intrinsic degrees of freedom discussed above. 



Indeed, the semiclassical quantization of the baryon current in Eq. (|2JJ) leads the appearance of 
the terms proportional to J 2 /M^ and I 2 /Mn, where is the nucleon mass. These terms represent 
coupling between the collective and intrinsic excitations. However, as discussed above these terms 
are of order l/N c and do not contribute at leading order in the l/N c expansion. This result allows 
one to find expectation values of the inclusive cross section using the superposition of the initial 
hedgehogs described by the collective variables A\ and A2 weighted by the corresponding Wigner 
D-matrix. In other words, the spin-isospin part of the nucleon wave function in the initial state with 
given quantum numbers m^m 1 is 

\m J ,m I > = J dA D]^ (A) \A > , (2.2) 

where [A > represents a hedgehog with particular orientation in spin-isopsin space, D 1 / 2 is the 
Wigner D-matrix describing spin-isospin coordinates of the nucleon, and the integral is taken over 
the space of the collective coordinates. 



Using Eq. one can find the inclusive (integrated over all mesons in the final state) nucleon- 



nucleon differential cross section at leading order in 1 /N c expansion, 

dat jffisM = [ dAldAl[D 'v ,( Al) iV /2 , w p* h( y i * 1 , <w> 

dQ. J m i' m i m ^ m 2 dQ ' 



where the dc mc (v,A\ ,A 2 ; 6, ty)/dQ is given in Eq. (2T). In Eq. (23), m[,m I l ,m2,m I 2 are the spin 
and isospin components of two nucleons along the direction n which can be taken as the beam axis. 
It is possible now to find a general form of the inclusive differential cross section at leading 



order in 1 /N c by integrating over the impact parameter space in Eq. ( [2T[ ) and SU (2) measure in 
Eq. (p3|). The resulting expression found in [P is 

do ( - mi '< m2 >^(v,0,(l>)/d£l = 
a o (v,e,0)+&/(v,e,0) (di -a 2 ) (?i-T 2 ) + c/O,0,0) (di ■n)(d 2 -n) (?i-T 2 ) , (2.4) 

where a' and x a are the spin and isospin Pauli matrices corresponding to the two initial nucleons, 
and the functions ciq, bj and cj encode the leading order behavior at large Af c . In obtaining the result 



in Eq. the following identity of the Wigner D-matrices was used, 

(Di hn y = (-l) m -"Di nh _ n . (2.5) 
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Equation (2.4) represents an inclusive differential cross section at leading order in 1 /N c expan- 
sion. However, more readily available is the data for total inelastic nucleon-nucleon cross section. 
To obtain the total cross section one should integrate the differential cross section over the whole 
solid angle. In doing so one obtains the following expression, 

aMmrtify) = Ao (v) +B/(v) (a, • a 2 ) (?, • ? 2 ) + C,(v) (a, - ft) (a 2 -n) ft • f 2 ) , (2.6) 

where 

Ao(v) = 1/2 J dQ.a (v;6,<l>), fl 7 (v) = 1/2 / dabr(v,e,^), C/(v) = 1/2 / dQ, Cl {v\d^). 

The factor of 1/2 in the above equations are due to the normalization of the baryon current. 

However, while formally integrating over the solid angle to obtain Eq. (E^) we did not consider 



that as discussed above the TDMFT treatment from which the differential cross section in Eq. (2.3) 



was derived strictly holds only in the semiclassical limit. As is well-known, Q19Q the semiclassical 
approximation breaks down for small forward angles for which 6 < 1 /Rpn, where R ~ N% is a size 
of the interaction potential which is of the order of the nucleon size, and ~ N c is the center-of- 
mass momentum. Thus, the forward angle at which the semiclassical approximation breaks down 
is of order of l/N c in the kinematic region of interest here. However, the total cross section in 



Eq. (2.6) includes both forward and backward angles for which the semiclassical approximation 
breaks down. 

However, as shown in details in [jl|] the contribution to the total cross section from the forward 
angles 6 < l/N c vanishes as l/N c in the large N c limit provided the scattering cross section is not 
anomalously peaked in a vanishingly small forward direction. The latter would require large elastic 
contribution. Thus, the total cross section given in Eq. (^^) is valid up to corrections of order N c . 

The key result of the above discussion is the form of total nucleon-nucleon cross section given 
in Eq. ( ^^ ) which is valid up to corrections of order 1 /N c . While the Skyrme model has been used 
in the derivation, the result is independent of the details of the model and are valid in the large N c 
limit. All the details of the dynamics for which the explicit TDMFT equations have to be used are 
in the functions A$(v), Bq(v) and Cb(v) which are of order one in the large N c limit. These functions 
can not be determined explicitly at present time. Nevertheless, the form of the total cross section in 



Eq. (2.6) does contain testable predictions. Note that it is not the most general form which can be 
obtained based on parity and time reversal invariance. For example, it does not contain such terms 
as A[ ft • f 2 ) and B ■ a 2 ). These terms do not appear at leading order in 1 /N c . This result 
does not depend on the details of the Skyrme but follows from the contracted spin-flavor symmetry 
SU (4) applied to the two-nucleon sector. 

This result can be tested against the existing experimental data for total nucleon-nucleon cross 
sections. This is done in the next section [p. 

3. Comparison with experimental data 

The experimental data exists for the total spin-independent and polarized proton-proton and 



proton-neutron cross sections [20]. The kinematic regime in which the result obtained above is 



expected to be valid corresponds to the center-of-mass momentum above 1 GeV . 
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Figure 1: Spin-averaged proton-proton and neutron-proton total cross section as a function of beam mo- 
mentum (Bugg et al., 1996). 



The total nucleon-nucleon cross section given in Eq. ( |2.6| ) contains the total nucleon-nucleon 
cross sections for isosinglet and isotriplet initial states, 



-(mi ,nt2 ,m 2 ) 



(v) 



<y(/=0) 

a( /=1 ) 



Using projection operators (1 — Z\ ■ T2) /4 and (3 — T1 • T 2 ) /4, one can extract a' /=0 ) and cj( /=1 ) 
cross sections. Then one obtains the total proton-proton, a pp = a^ I=1 \ and neutron-proton cross 

Thus, at leading order in l/N c one has the following expres- 



sections, a np 
sions, 



r(/=l) 



+ 



t {pp) 



(1=0)) 



A Q (cti • ct 2 ) + Q (cti • n) (ct 2 • n) 
Aq-Bj (di -d 2 )-C/ (3i ■n)(a 2 -n) . 



(3.1) 



Recall that in the above equations H represents the beam axis. Using Eq. (3.1 ) and averaging over 
the spin-polarization of the beam and target nucleons, one can obtain spin-averaged total cross 
sections. Thus, at leading order we have the following relation, 



^ = ^(1 + ^(1/^)), 



(3.2) 



where Ob's are the spin-averaged total cross sections. This prediction follows from the large- 
A^ c analysis and cannot be obtained simply from isospin invariance. This large-A^ c result is well 
satisfied by the experimental data shown in Fig. [1} 

The total cross sections for the case when beam and target nucleons are transversely polarized 
relative to the beam direction can also be obtained. Two configurations are possible, \\ and tl- 
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Figure 2: Ao^ for neutron-proton (Fonteine et al., 1991) and proton-proton (Ditzler et al., 1983; Lesikar, J. 
D. 1981) scattering, and Aol for proton-proton (Auer et al., 1978), and neutron-proton (Sharov et al., 2008) 
scattering as functions of beam momentum. 



These can be combined into an observable, Aaj = — (ct(Tt) ~~ ^(tl)) referred to as delta sigma 
transverse. 

Analogously, for the the longitudinally polarized beam and target nucleons one can extract an 
observable, Acl = — (cr(=t) — a(<^)), the delta sigma longitudinal. 



The large-A^ analysis, Eq. ( p.l| ), predicts the following results for these observables, 



Ac 
Act. 



(pp) 

T 

(PP) 



-Ao!t np) (l + 0(1/N c )), 
-Aol" p) (\ + 0(\/N c )). 



(3-3) 



Experimental data for these observables is shown in Figs. [2(a)| and [2(b)[ One may conclude 



from Figs. |2(a)| and [2(b)| that the large N c results given in Eq. ( p.3| ) are not satisfied by data 



However, the results are valid within corrections of order N c . Indeed, according to Eq. ( |2.6[ ) both 
Co and Aot and Actl are of the same order in 1 /N c . Experimentally however Act^ and Actl are 
much smaller then Ob. The latter are about 40mb while the former are consistent with zero. The 
suppression of Act^ and Actl are for reasons not predicted by large-A^ analysis. Nevertheless, 
qualitatively the predictions are valid since ACT7- and Aol are small for both pp and np scattering. 
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